Shape transition and oblate-prolate coexistence in N = Z fpg-shell nuclei 
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Nuclear shape transition and oblate-prolate coexistence in TV = Z nuclei are investigated within 
the configuration space (2p 3 / 2 , I/5/2, %Pi/2, an d 159/2). We perform shell model calculations for 
60 Zn, 64 Ge, and 68 Se and constrained Hartree-Fock (CHF) calculations for 60 Zn, 64 Ge, 68 Se, and 
72 Kr, employing an effective pairing plus quadrupole residual interaction with monopole interactions. 
The shell model calculations reproduce well the experimental energy levels of these nuclei. From the 
analysis of potential energy surface in the CHF calculations, we found shape transition from prolate 
to oblate deformation in these N = Z nuclei and oblate-prolate coexistence at 68 Se. The ground 
state of 68 Se has oblate shape, while the shape of 60 Zn and 64 Ge are prolate. It is shown that the 
isovector matrix elements between / 5 / 2 and P1/2 orbits cause the oblate deformation for 68 Se, and 
four-particle four-hole (4p — ih) excitations are important for the oblate configuration. 

PACS numbers: 21.60.Cs, 21.60.Jz, 21.60.-n,21.10.-k 
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The proton-rich nuclei with masses between 56 and 80 
are well known to change rapidly their collective prop- 
erties with proton and neutron numbers. The relevant 
Nilsson diagram pj shows a large shell gap for N, Z =34 
at prolate and oblate deformations, for N, Z =36 at large 
oblate deformation, and for N, Z =38 at large prolate 
deformation. Particular interests are in shape transition 
and oblate-prolate coexistence of heavier N = Z nuclei. 
These nuclei lie in the transitional region from the spher- 
ical shape (e.g. 56 Ni 0) with coexisting prolate shape 
Q to the prolate deformation (e.g 
cleus 64 Ge p| is known to be an N 
unstable nucleus manifesting a 7-soft structure from the 
theoretical calculations based on the mean-field approx- 
imation. The early studies of 69 Se |(J have found indi- 
cations of oblate shapes. Long standing predictions of 
a stable oblate deformation were recently confirmed by 
the observation of an oblate ground state band in 68 Se 
. Determinations of shape were inferred indirectly from 
the study of rotational bands, while direct quadrupole 
measurements are difficult for short-lived states. In their 
analyses, it was suggested that the oblate configuration 
coexists with a prolate rotational band which is so-called 
oblate-prolate coexistence. Previous work Q for 72 Kr 
suggests that the ground state has an oblate shape and 
there is an oblate — ► prolate shape transition at low spins. 
In the total Routhian surface (TRS) calculations 0, the 
oblate minimum can be understood in terms of the large 
gap for oblate deformation at /V = Z =36. Projected 
shell model a deformed selfconsistent Hartree-Fock 
(HF)+RPA approach and Skyrme Hartree-Fock- 

Bogolyubov (HFB) calculations [12j f° r proton-rich nu- 
clei in the A =70-80 mass region were performed. Recent 
investigations for 68 Se have been carried out by the pro- 
ected shell model |l3j, by the excited VAMPIR method 
141 , and the self-consistent collective coordinate method 

HI- 

The spherical shell model approaches could be more 



appropriate for describing various aspect of nuclear struc- 
ture. However, the shell model calculations by diagonal- 
ization in the l/ 7/2 , 2p 3/2 , l/ 5/2 , 2p 1/2 , and lg g/2 for nu- 
clei with N, Z = 30 — 36 are hopeless at present because 
of huge dimension of configuration space. So we need 
to restrict the model space to the 2p 3 / 2 , 1/5/21 2pi/ 2 , and 
l#9/2 orbits (henceforth called fpg-sheW). The dimension 
of the configuration space is still huge. For instance, max- 
imum dimension for 68 Se is 0.165 billion. Recently, an 
extended pairing plus quadrupole-quadrupolc (P + QQ) 
force has been applied to the fpg-sheW nuclei, and shown 
to be useful [la, [jjl • This interaction works remarkably 
well. In this paper, we study the shape transition and 
oblate-prolate coexistence in the N = Z fpg-sheW nu- 
clei using a large scale shell model and the constrained 
Hartree-Fock (CHF) calculations. The shell model cal- 
culations with ~ 10 8 dimension can be carried out by 
recently developed shell model code jig. The observa- 
tion associated with oblate shape is rare compared with 
the prolate shape well established experimentally. The 
reason for the suppression of oblate deformation lies in 
the higher order effects both in liquid drop terms and 
residual interactions which make most nuclei favor pro- 
late shapes. Very strong oblate-driving effects seem to 
be necessary for overcoming this prolate tendency. As 
a candidate for the oblate-driving force, we present the 
isovector (T = 1) interaction with matrix element be- 
tween / 5 / 2 andp!/ 2 orbits. 

We start from the extended P 
with the monopole interactions V m 
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where e a is a single-particle energy, Pooik is the T = 
1, J = pair operator, and Q 2 m {O^m) is the isoscalar 
quadrupole (octupole) operator. Each term includes p-n 
components which play important roles in the N = Z 
nuclei, due to an isospin-invariance. 

The shell model calculations pH have been per- 
formed in a shell model space, which assumes a closed 
56 Ni core. The neutron single-particle energies of 2p 3 / 2 , 
I/5/2, 2p 1/2 , and lg 9/2 in this fpg-shell region can be 
read from the low-lying states of 57 Ni. We used the mea- 



sured values £ P 3/2 = 0.0, 2/5/2 = 0.77, £ p i/2 = 1.11, and 
e g9/2 = 3.70 in MeV, in our previous paper 0. How- 
ever, we have recently obtained a better single particle 
energy £ m/? = 2.50 MeV which fits odd-mass Ge isotopes 
in Refs. |2Ct l2l| . This value is consistent with the report 
in [2^ • Since the above Hamiltonian is assumed to be an 
isospin-invariant, the proton single-particle energies are 
taken as the same values as the neutron single-particle 
energies. We searched force strengths of the extended 
P + QQ interaction using £ g9 / 2 — 2.50 MeV. The ob- 
tained values are as follows: 



go = 0.270(64/A), X 2 = 0.250(64/A) 5/3 /6 4 , 
X3 = 0.05(64/A) 2 /6 6 , fc° = 1.44(64/A), 



(4) 



where go, X2, X3i an d ^° are the J = pairing, the 
QQ, the octupole, and J-independent isoscalar force 
strengths, respectively. Here b is the harmonic-oscillator 
range parameter. The monopole shifts are 

V m (p3/2,f 5 /2;T = 1) = V m (p 3/2 ,p 1/2 ;T = 1) = -0.3, 

Kn(/5/2,Pl/ 2 ;T=l) = -0.4, 

V m (g 9 /2,g 9 / 2 ;T = 1) = -0.2, 

V m (g 9 /2, 59/2; T = 0) = -0.1, in MeV (5) 

All these force strengths have been phenomenologi- 
cally adjusted so as to fit many energy levels including 
high-spin levels for 60,62,64,66,68 Zllj 64,66,68,70 Ge and also 

65 ' 67 Ge ji3,|2l|. The force strengths determined in this 
way describe quite well various properties of these nuclei. 
In particular, the monopole shifts are important for de- 
scribing precise positions of the high-spin (J > 8) states 
for 66 > 68 Ge. 

In Fig. 1, energy spectra calculated with the above 
force strengths are compared with experimental data in 
the N = Z fpg-shell nuclei, 60 Zn, 64 Ge, and 68 Se. In 
addition to the ground state band, the second positive- 
parity band beginning from J" = 2 2 state and negative- 
parity band beginning from J 7 ' = 3j~ state are shown 
in Fig. 1. The calculated energy levels are in good 
agreement with the experimental energy levels, except 
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FIG. 1: Experimental and calculated energy levels of 
64 Ge, and 68 Se. 



J Zn, 



for r = 4+ and r 



Q 2 levels for 64 Ge. In particular, 



the energy levels of 68 Se are well reproduced. The exper- 
imental level sequence in the ground-state band shows a 
quadrupole deformation. 

In order to study the quadrupole deformation of these 
nuclei, we examine the spectroscopic quadrupole mo- 
ments and E2 transitions. We adopted the effective 
charge e p = 1.5e for proton and e„ = 0.5e for neutron. 
Figure 2 shows the calculated spectroscopic quadrupole 
moments for the first and second excited 2 + states for 
60 Zn, 64 Ge, and 68 Se. The quadrupole deformation was 
estimated from the spectroscopic quadrupole moments 
for the first and second excited 2 + states on the assump- 
tion of axially symmetric deformation. For 60 Zn and 
64 Ge, the deformation estimated from the quadrupole 
moments are prolate with (3 ~ 0.2 for the first excited 
2 + state (2^) and oblate with /3 ~ —0.2 for the second 
excited 2 + state (2 2 ). The nucleus 68 Se favors an oblate 
deformation with j3 ~ —0.2 since the value of intrinsic 
quadrupole moment for the 2^~ state becomes negative. 
On the other hand, the 2 2 state has prolate deformation 
with (3 ~ 0.2. The results indicate shape transition in 
these nuclei. We examined the contributions of interac- 
tion matrix elements to the quadrupole moments. We 
found that the T = 1 matrix elements between / 5 / 2 and 
Pi/ 2 orbits strongly affect the quadrupole moments. The 
T = 1 monopole interaction is particularly important for 
the inversion of signs of the quadrupole moments in 68 Se. 
As seen in Fig. (2), the inversion of the quadrupole mo- 
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FIG. 2: Spectroscopic quadrupole moments in Zn, Ge, 
and 68 Se. The solid squares and circles are quadrupole mo- 
ments of the first and second excited 2 + states, respectively, 
obtained by the shell model calculations. The open squares 
and circles are those obtained by ignoring the T = 1 monopole 
matrix elements. 



lated by Petrovici et al. 44]. This can be attributed to 
neglection of excitations from the /V/2 orbit because of 
computational limitation. 

Let us next examine the nuclear shapes including tri- 
axiality of 60 Zn, 64 Ge, 68 Se, and 72 Kr by an alternative 
approach, the CHF method which is carried out by 
the following quadratic way: 



H' 



H 



) - qrf + (3{{J X ) - 3x )\ (6) 



where J x are the x-component of the angular momen- 
tum operator. The q^s are the constant parameters 

9o = v/db? cos 7> <7±2 = \/T^9 sin 7. and Q±i = 0, where 
q is the isoscalar intrinsic quadrupole moment and 7 is 
the triaxial angle. We set j x — J(J + 1) with J being 
the total angular momentum of the state. The param- 
eters, a and /3, are taken so as to achieve convergence 
of the iteration of the gradient method. Then, poten- 
tial energy surface (PES) is denned as the expectation 
value (H) with respect to the CHF state for a given q 
and 7. Figure 3 shows the contour plots of the PES in 
the q-7 plane for 60 Zn, 64 Ge, 68 Se, and 72 Kr. We can see 
remarkable shape changes of these nuclei. The PES mini- 
mum exhibits prolate deformation for 60 Zn and 64 Ge, and 
oblate one for 68 Se, and triaxial one for 72 Kr. The char- 
acteristic features of the PES's are the 7-softness for 64 Ge 



ments occurs in 68 Se, when we remove only the T = 1 
monopole matrix elements v = V m (f^/2iPi/2\T = 1) = 
—0.4 MeV from all of the interaction matrix elements. 
Namely, the shape of ground state in 68 Se becomes pro- 
late in this case. Thus, we can say from our calculations 
that the inverse signs of the quadrupole moments are at- 
tributed to the T = 1 interaction with matrix element 
between / 5 / 2 and px/2 orbits. 

The quadrupole deformation of Ge estimated from 
the calculated B(E2;2+ -> 0^) is (3 ~ 0.2 16]. This 
value is in agreement with the predictions (3 ~0.22 by 
Moller and Nix [23|, and (3 - 0.22 by Ennis et al. Of. 
The calculated value of B(E2; l\ -> 2\)/B(E2; 2\ o\) 
is very large, which indicates the 7-softness according to 
the Davydov model [24j|. This is in agreement with the 
7-softness or triaxiality estimated from the experimental 
data and the other theoretical models. It was demon- 
strated p^ ] that the proton-neutron part Q p Q n of the 
quadrupole-quadrupole interaction is important for the 
7-softness or triaxiality. For 68 Se, the quadrupole defor- 
mation estimated from the calculated B(E2; 2f — * Of) 
is (3 ~ 0.21. The deformations estimated from the 
quadrupole moment and B(E2) are smaller than (3 ~ 
0.27 estimated by Fischer et al. 0, P ~0. 30 obtained by 
Jenkins et al. [25j for 69 Se and the value (3 ~ 0.33 calcu- 
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FIG. 3: Contour plots of PES's on q-7 plane in the CHF 
calculations for 60 Zn, 64 Ge, 68 Se, and 72 Kr. The solid circles 
denote the PES minima, and the solid triangle the second 
minimum. 
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and oblate-prolate coexistence for 68 Se,which are consis- 
tent with the previous discussions [E0- Moreover, we 
calculated an angular-momentum projected PES which is 
obtained by calculating the expectation value (H) with 
respect to the angular-momentum projected CHF wave 
function for a given q and 7. The angular-momentum 
projection method improves considerably the ground- 
state energies. The angular-momentum projected PES's 
show a similar potential energy surface to the CHF. 

The PES's as a function of q along the axially symmet- 
ric line are shown for 60 Zn, 64 Ge, 68 Se, and 72 Kr in Fig. 
4. This figure shows the aspect of shape coexistence and 
shape changes in these nuclei, from different angle. Fig- 
ure 4 also shows that if the T = 1 monopole interactions 
are removed, the minima of the PES's in 68 Se and 72 Kr 
shift to the prolate side, while it remains unchanged in 
60 Zn, 64 Ge. Thus the T = 1 monopole interactions with 
the matrix element V m {f5/2,Pi/2\ T = 1) are confirmed 
to be the oblate-driving force for 68 Se and 72 Kr. The 
shape of the ground state becomes oblate for 68 Se where 
the Fermi energy for neutron lies between the / 5 / 2 and 
P1/2 orbits. The CHF calculation shows that the shape 
of 72 Kr is triaxial, while it is oblate in the HFB result 
[13 and the TRS results @. Note that the PES of 72 Kr- 
is flat toward oblate shape and shows gamma softness, 
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FIG. 4: PES's as a function of q along the axially symmetric 
line. The solid curves are the PES's in the full shell model 
matrix elements and the dotted curves are those obtained by 
ignoring the T = 1 monopole matrix element between / 5 / 2 
and p 3/2 . 
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FIG. 5: Spectroscopic quadrupole moment as a function of 
t. The solid squares and circles are for the first and second 
excited 2 + states, respectively. The open squares and circles 
are those obtained by ignoring the T = 1 monopole matrix 
elements. 



as seen in Fig. 3. 

We next analyze how excitations from the orbits 
(P3/2>/5/2) to (pi/2,59/2) contribute to the oblate de- 
formation. We calculate the quadrupole moments by 
enlarging configuration space in the truncation scheme 
s <t(P3/2/5/2) j4_56 " s (Pi/259/2) s , where t is the maxi- 
mum number of particles allowed to be excited. Figure 5 
shows the spectroscopic quadrupole moments as a func- 
tion of t for 68 Se. The spectroscopic quadrupole mo- 
ment is negative for 2^~ for l\ when t < 3, but suddenly 
their signs change at t=A. This means that the four- 
particle four- hole (4p — Ah) excitations from (P3/2/5/2) 
to (pi/239/2) are very important for the oblate configu- 
rations. Such a change does not occur when the T = 1 
monopole interactions are removed, as shown by the dot- 
ted lines in Fig. 4. Table I shows the occupation numbers 
of the 2p 3 / 2 , I/5/2) 2pi/2, and 1^9/2 orbits for the first 
and second excited 2 + states of 60 Zn and 68 Se in the 
cases of the T = 1 monopole force strengths v = —0.4 
MeV and 0.0 MeV. When we add the T = 1 monopole 
force, the occupation number of p\/2 for the first excited 
2 + state for 68 Se increases remarkably. Thus the T = 1 
monopole force contributes to Ap — Ah excitations from 
(P3/2/5/2) to P1/2, which is especially important for the 
oblate configuration. 

Our conclusion is different from that of the VAMPIR 
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TABLE I: Occupation numbers in the first and second excited 
2+ states for 60 Zn and 68 Se. 



60 Zn 


2f 
l 

P3/2 /5/2 Pl/2 39/2 


2+ 

2 

P3/2 /s/2 Pl/2 39/2 


v=-0A 
v= 0.0 

68 Se 


0.98 0.51 0.49 0.02 
1.06 0.46 0.46 0.02 

P3/2 /5/2 Pl/2 59/2 


1.08 0.59 0.24 0.09 
1.11 0.56 0.24 0.09 

P3/2 /5/2 Pl/2 39/2 


v=-0A 
v= 0.0 


2.31 2.38 1.20 0.11 
2.65 2.78 0.42 0.15 


2.37 2.81 0.69 0.13 
2.70 2.33 0.80 0.15 



calculations [l4|, where the oblate-driving force for 68 Se 
is in monopole shifts for T = matrix elements between 
the g 9 /2 and / 5/2 (/7/2) orbits. We examined this question 
by changing monopole shifts and single-particle energy 
e g g /2 . The investigations of the Q moments show that the 
monopole shift V m (/s/2, 99/2! T = 0) hardly contributes 
to the oblate deformation under the single-particle ener- 
gies adopted in this paper. If we lower e g g/2 below e p i/2 
as done in Rcf. , many nucleons occupy the gg /2 orbit 
and V m (f 5 / 2 , gg/2'1 T = 0) contributes to the oblate defor- 
mation. The nucleon occupation numbers of the P1/2 and 
gg j2 orbits and phases of respective components of wave- 
functions affects the sign of the Q moments. It seems 
that the VAMPIR conclusion depends on the inversion 
of £p 1/2 and e g g/2- 

Finally, we discuss shortly standard P + QQ force 
model calculations neglecting the octupole force and the 
monopole shifts. We set the force parameters such as so 
as to reproduce the energy levels of ground-state band for 
68 Se. Note that the quadrupole force strength is larger 
than that in Eq. (4). Then, we obtain the positive spec- 



troscopic quadrupole moment for the 2f state and neg- 
ative one for the 2% state. Similarly in the CHF calcu- 
lations, the first PES minimum is in the oblate side and 
the second minimum is in the prolate side. However, the 
standard P + QQ force model cannot consistently repro- 
duce the energy levels of the side band for 68 Se and the 
ground-state bands for 60 Zn and 64 Ge. 

In conclusion, we investigated nuclear shape transition 
and oblate-prolate coexistence in N = Z /py-shell nu- 
clei. We analyzed the quadrupole moments in the shell 
model and the CHF calculations. The analysis shows 
that the T = 1 monopole interactions with the matrix ele- 
ment between / 5 / 2 and P1/2 orbits play an important role 
in producing the oblate shape for 68 Se and 72 Kr, while 
they do not play the same role for 60 Zn and 64 Ge. Note 
that this conclusion is different from that of the VAM- 
PIR calculations by Petrovici et al., where the origin of 
oblate-driving force is T — monopole matrix elements 
between gg/ 2 and j '5/2(77/2) orbits. However, in our cal- 
culation we could not find any large influence from these 
T = monopole matrix elements on the oblate deforma- 
tion. The shape transition and oblate-prolate coexistence 
in the CHF calculations are consistent with those of the 
Skyrme HFB calculations The present calculations 
show that the oblate deformation originates in the Ap—Ah 
excitations from (P3/2/5/2) to (^1/239/2)- We can expect 
that the T = 1 monopole interactions between / 5 / 2 and 
P1/2 orbits are also important for the oblate or triaxial 
deformation in N 7^ Z nuclei where the Fermi surface for 
neutron lies between / 5 / 2 and P1/2 in fpg shell. Further 
investigations are in progress. 
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